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Abstract 

The highest energy cosmic rays, above the Greisen-Zatsepin-Kuzmin cut-off 
of cosmic ray spectrum, may be produced in decays of superheavy long-living 
X-particles. We conjecture that these particles may be produced naturally in 
the early Universe from vacuum fluctuations during inflation and may consti- 
tute a considerable fraction of Cold Dark Matter. We predict a new cut-off 
in the UHE cosmic ray spectrum E cut - ff < mi n fl a t on ~ 10 13 GeV, the exact 
position of the cut-off and the shape of the cosmic ray spectrum beyond the 
GZK cut-off being determined by the QCD quark/gluon fragmentation. The 

Pierre Auger Project installation might discover this phenomenon. 
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According to the Greisen-Zatsepin-Kuzmin [|TJ (GZK) bound, the Ultra High Energy 
(UHE) cosmic rays produced in any known candidate extra galactic source should have an 
exponential cut-off at energies E ~ 5 x 10 10 GeV. On the other hand, the number of observed 
cosmic rays events beyond the cut-off is growing and leads to a mounting paradox within 
standard frameworks of cosmological and particle physics models. 

A wide variety of possible solutions were suggested. Resolution could be due to exotic 
particle which may be produced at cosmological distances were suitable conventional ac- 
celerators are found, be transmitted evading GZK bound, and yet which interact in the 
atmosphere like a hadron. A particle with correct properties was found in a class of su- 
persymmetric theories ||. Alternatively, high energy cosmic rays may have been produced 
locally. One possibility is connected to the events of destruction of (topological) defects [|J, 
while another one to decays of primordial heavy long-living particles |5|||. The candidate 
particle must obviously obey constraints on the mass, density and lifetime. 

In order to produce cosmic rays of energies E > 10 11 GeV, the mass of X-particles has 
to be very large, mx ^ 10 13 GeV [§,§J. The lifetime, tx, cannot be much smaller than the 
age of the Universe, r > 10 10 yr. With this smallest value of the lifetime, the observed flux 
of UHE cosmic rays will be reproduced with rather low density of V-particles, Qx ~ 10~ 12 , 
where Qx = "^x^x/Pcrit, nx is the number density of X-particles and p cr i t is the critical 
density. On the other hand, X-particles must not overdose the Universe, Qx ~ 1- With 
Qx ~ 1 , the X-particles may play the role of cold dark matter and the observed flux of UHE 
cosmic rays can be matched if Tx ~ 10 22 yr. The allowed windows are quite wide ||, but 
on exotic side, which may rise problems. 

The problem of the particle physics mechanism responsible for a long but finite lifetime 
of very heavy particles can be solved in several ways. For example, otherwise conserved 
quantum number carried by the X-particle may be broken very weakly due to instanton 
transitions |J, or quantum gravity (wormhole) effects |J. If instantons are responsible for 
V-particle decays, the lifetime is roughly estimated as Tx ~ Tn x • exp(47r/ax), where ax 
is the coupling constant of the relevant (spontaneously broken) gauge symmetry. Lifetime 
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will fit the window if the coupling constant (at the scale mx) is ax ~ 0.1 0. 

X-particles can be produced in the right amount by usual collision and decay processes 
if the reheating temperature after inflation never exceeded mx, but the temperature should 
be in the range 10 11 ^ T r < 10 15 GeV, depending on mx, This is a rather high value 

of reheating temperature and will lead to the gravitino problem in generic supersymmetric 
models [[?]]. 

In the present paper we investigate another process of X-particle creation, namely the 
direct production from vacuum fluctuations during inflation. 

Any viable modern cosmological model invokes the hypothesis of inflation ||. During 
inflation the Universe expands exponentially, which solves the horizon and flatness problems 
of the standard Big-Bang cosmology. Inflation is generally assumed to be driven by the 
special scalar field <fi known as the inflaton. Fluctuations generated at inflationary stage can 
have strength and the power spectrum suitable for generation of the large scale structure. 
This fixes the range of parameters in the inflaton potential. For example, the mass of the 
inflaton field has to be ~ 10 13 GeV. During inflation, the inflaton field slowly rolls down 
towards the minimum of its potential. Inflation ends when the potential energy associated 
with the inflaton field became smaller than the kinetic energy. At that time all the energy 
of the Universe is contained entirely in the form of coherent oscillations of the inflaton field 
around the minimum of its potential. It is possible that a significant fraction of this energy 
is released to other Boson species after only a dozen or so inflaton field oscillations, in the 
regime of a broad parametric resonance ||. This process was studied in details [pT0| , pLlJ1 . Even 
particles with masses of order of magnitude larger than the inflaton mass can be produced 
quite abundantly. Applying these results to the case of our interest here, we find that stable 
very heavy particles, mx ^ m^, generally will be produced in excess and will overdose the 
Universe. 

However, if the parametric resonance is ineffective for some reason, and we estimate 
particle number density after inflation at the level of initial conditions used in Refs. [|10|| , we 
find that fix might prove to be of about the right magnitude. This level is saturated by the 



fundamental process of particle creation during inflation from vacuum fluctuations and it is 
the same process which generated primordial large scale density fluctuations. Parametric 
resonance for X particles is turned-off if X is either a fermion field or its coupling to inflaton 



is small, g 2 < 10 4 (m x M/,) ("W 'M Pl ) 2 |0|. 



At some early epoch the metric of the Universe is conformally flat to a high accuracy, 
ds 2 = a(n) 2 (dn 2 — dx 2 ). We normalize the scale factor by the condition a(0) = 1 at the end 
of inflation. Number density of particles created in time varying cosmological background 
can be written as 

nx = ^J m 2 k 2 dk , (1) 

where /3& are the Bogoliubov coefficients which relate "in" and "out" mode functions, and k 
is the comoving momentum. Massles conformally coupled particles (for scalars this means 
that £ = 1/6 in the direct coupling to the curvature) are not created. For massive particles 
conformal invariance is broken. Therefore, for the power low (e.g., matter or radiation 
dominated) period of expansion of the Universe, one expects on dimensional grounds, nx oc 
m\ j a? at late times. Indeed, it was found in Ref . [|EJ] 

n x ~ 5.3 x 10~ 4 m x (m x t)' 3/2 , (2) 

for the radiation dominated Universe, and nx oc m x (m x t)~ 3q for a(t) oc t q . Note that all 
particle creation occur in the region rat = qm/H < 1. When mt <C 1, the number density 
of created particles remains on the constant level nx = m x /24ir 2 independent of q []12| |. At 
qm/H 3> 1 particle creation is negligible. Here H is the Hubble constant, H = a/a. 

For the radiation dominated Universe one finds, Qx ~ (m^- / Mpj) y/ mxt e , where t e is 
time of equal densities of radiation and matter in fi = 1 Universe. This gives fix ~ m^ 2 , 
where mg = m^/lO 9 GeV. Stable particles with mx ~ 10 9 GeV will overdose the Universe 
even if they were created from the vacuum during regular Friedmann radiation dominated 
stage of the evolution. It is possible to separate vacuum creation from creation in collisions 
in plasma since X-particles may be effectively sterile. 



This restriction can be overcomed if evolution of the Universe, as it is believed, was more 
complicated than simple radiation dominated expansion from singularity. Hubble constant 
may have never exceeded mx , which is the case of inflation, H(0) ~ m^. Moreover, compared 
to the case considered above, density of X-particles created during inflation is additionally 
diluted by late entropy release in reheating after inflation. 

Particle creation from vacuum fluctuations during inflation (or in de Sitter space) was 
extensively studied fl^JT^I . Characteristic quantity which is usually cited, the variance of 
the field (X 2 ), is defined by an expression similar to Eq. ([!]). In the typical case ~ — /3& 
the difference is given by the factor 2 sin 2 (u^r/) /uk in the integrand, where uj\ = k 2 + a 2 m 2 x . 
If nix ~ H(0) m m^, one has on dimensional grounds rix = Cm^/Zir 2 ^ where the coef- 
ficient C is expected to be somewhat smaller than unity. Both Fermions and Bosons are 
prodused by this mechanism, exact numerical value of C being dependent on spin-statistics. 
In general, C is the function of the ratio H(0)/mx, the function of self-coupling of X and 
the coupling £, depends on details of the transition between inflationary and matter (or radi- 
ation) dominated phases, etc. For example, for the scalar Bose-field, (X 2 ) = 3H(0) /8ir 2 Tn x 
if m x < H(0) 0,0. For massless self-interacting field (X 2 ) « 0.132/7 (O) 2 /^ |§. C is 
expected to decrease exponentially when m x > m^. Particle creation in the case of Hubble 
dependent effective mass, mx(t) oc H(t), was considered in Ref. [|l(| . 



Let us estimate the today's number density of X-particles. We consider massive inflaton, 
V{4>) = m 2 .0 2 /2. In this case inflation is followed by the matter domination stage. If there 
are light Bosons in a theory, < m^, even relatively weakly coupled to the inflaton, 
g 2 > lC^m^/Mp! ~ 10~ 8 , this matter domination stage will not last long: inflaton will decay 
via parametric resonance and the radiation domination follows. This happens typically when 
the energy density in inflaton oscillations is redshifted by a factor r w compared to 
a value m 2 h M^, l Matter is still far from being in the thermal equilibrium, but it is 

still convenient to characterize this radiation dominated stage by an equivalent temperature, 



T* ~ r 1 / 4 ym^Mpi. At this moment the ratio of energy density in X-particles to the total 
energy density retains its value reached at the end of inflation, p x / Pr ~ C m^mx /27r 2 Mpj. 



Later on this ratio grows as oc T/T* and reaches unity at T = T eq , where 

CW4/ m x3/2 
T eq = TT- ™* • ( 3 ) 
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Using relation T eq = 5.6fi x /i 2 eV we find that 1(T 12 < fi x £ 1 if 

1(T 23 < Cr^mx/m* < 10" 11 . (4) 

For rax ~ (a few) • this condition can be easily satisfied since the coefficient C is 
exponentially small. This condition may be satisfied even for rax ~ ra<t> since the coefficient 
r 1//4 (or equivalent reheating temperature) can be small too. 

Our hypothesis has unique observational consequences. If UHE cosmic rays are indeed 
due to decay of superheavy particles which were produced from vacuum fluctuations during 
inflation, there has to be a new sharp cut-off in the cosmic ray spectrum at energy somewhat 
smaller rax- Since the number density rix depends exponentially upon rax/ra^, the position 
of this cut-off might be well predicted and has to be near -E cu t-ofr < ra^ w 10 13 GeV, the 
very shape of the cosmic ray spectrum beyond the GZK cut-off being of quite generic form 
following from the QCD quark/gluon fragmentation. The Pierre Auger Project installation 
| I7J might prove to be able to discover this fundamental phenomenon. 

We conclude, observation of Ultra High Energy cosmic rays can probe the spectrum of 
elementary particles in its superheavy range and can be an additional opportunity (alongside 
with fluctuations in cosmic microwave background) to study the earliest epoch of the Uni- 
verse evolution, starting from amplification of vacuum fluctuations during inflation through 
fine details of gravitational interaction and down to physics of reheating. 

When our paper was at the very end of completion we became aware of the quite recent 
paper by Chung, Kolb and Riotto QTSJ where similar problems of superheavy dark matter 
creation were considered. 
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